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Phonon modes and their dispersion relations in ultrathin homogenous dielectric membranes are 
' calculated using elasticity theory. The approach differs from the previous ones by a rigorous account 

f> | of the effect of the film surfaces on the modes with different polarizations. We compute the heat 

■ capacity of membranes and the heat conductivity of narrow bridges cut out of such membranes, 

in a temperature range where the dimensions have a strong influence on the results. In the high 
temperature regime we recover the three-dimensional bulk results. However, in the low temperature 
CI limit the heat capacity, CV, is proportional with T (temperature), while the heat conductivity, k, 

of narrow bridges is proportional to T 3 / 2 , leading to a thermal cut-off frequency f c = k/Cv oc T 1//2 . 



c3 



-a 
c 

o 
o 



> 

o 
o 

-I— > 



i 

o 
o 



X 



INTRODUCTION 

More and more precise measurements are needed in 
science (e.g. in astrophysics) and therefore very sensitive 
detectors are required. However, even if all the techno- 
logical difficulties are removed and the detector compo- 
nents would be perfect, the thermal fluctuations and the 
discreteness of particle and energy fluxes through the de- 
tector would still limit the precision of the measurements. 
To achieve the desired sensitivity, the detectors have to 
work at sub-kelvin temperatures (to minimize the fluc- 
tuations) and have linear dimensions in the micrometer 
scale or even smaller (for small heat capacity and for the 
possibility to assamble high resolution detector arrays for 
space applications). For such ranges of temperatures and 
dimensions, finite-size effects come into play. 

Good thermal insulation of the detector is essential 
to attain the working temperature. For this purpose, 
placing the detector on top of a free-standing, amorphous 
and dielectric silicon-nitride (SiN^) membrane appears 
to be a very convenient design (see for example Refs. 
0, S 0) Thc thickness of such a membrane is of the 
order of hundreds of nm, which is very small as compared 
to its other dimensions. The membrane may also be cut, 
so that the detector lies on its central and wider part, 
which is connected to the bulk material by narrow bridges 
(See Fig.Hl. 

Thermal characteristics of SiN^ membranes and 
bridges have been investigated experimentally by sev- 
eral authors (as examples, see Q, 0, H, @)- SiN^ is an 
amorphous dielectric and, if we neglect here the dynamic 
defects (see review Q and references therein), its low 
temperature thermal properties are determined only by 
phonons. At temperatures below a few degrees Kelvin, 
the bulk material (assumed to be isotropic) is well de- 
scribed by the Debye model, with typical values for the 
longitudinal and transversal sound velocities Ct = 6200 
m/s and c; = 10300 m/s, respectively 0. Hence, the 
heat capacity can be expressed as Cy D = KiT 3 , where T 






FIG. 1: (a) A schematic drawing of a thin free standing mem- 
brane of the type discussed in this paper. The thickness, b, is 
typically a few hundreds of nm. The other dimensions, h and 
h, are of the order of 100 /im. (b) Such a membrane (cen- 
tral part) suspended to the bulk by four bridges. The bridge 
width, w, ranges from few fim to tens of jim. 



is the temperature and K\ is a constant. If the mean free 
path of the phonons, £, is assumed to vary as I oc A~ s , 
where A is the phonon wavelength, the heat conductivity 
has the expression k 3D — KzT^ 8 , where again K2JS a 
material dependent constant (see for example Ref. @). 

Note that for the above given values of sound veloc- 
ity, and at temperatures below 1 K, the bulk dominant 
phonon wavelength, Adorn ~ he/ (1.6]cbT) (see Chap. 
8), is of several hundreds of nanometers. Thus it is com- 
parable with the smallest dimensions of the structures 
measured in Refs. Typically in these exper- 

iments, the measurements made at the lower end of the 
temperature interval showed qualitative differences from 
the measurements at the upper end of the interval, due 
to finite-size effects. 

In this paper we will improve the model of Refs. 
[ToL Hl| by calculating rigorously the phonon modes and 
their dispersion relations in ultrathin membranes and 
narrow bridges. We shall show the thermodynamical 
consequences of the calculations. Such calculations ap- 
pear to be well known to those working in the field of 
elasticity theory (see for example Refs. 0, El), but 
apparently are not so familiar to the condensed matter 
community. Therefore, to make the paper more read- 



able, we shall give some details of our calculations. To 
describe these effects, a simplified model was recently 
used 0, 0] . Namely, it was assumed that the phonon 
modes are quantized along the direction perpendicular to 
the membrane surface. In this way, the three-dimensional 
(3D) (quasi)continuous phonon spectrum splits into what 
one could call two-dimensional (2D) branches ^(J. At 
low enough temperatures, only the lowest branches are 
populated by phonons and the membrane is described 
as a 2D isotropic phonon gas. In this situation, if the 
mean free path of the phonons varies as I cx A _s , the 
heat capacity has the form Cy D — K3T 2 and the ther- 
mal conductivity is k 2D = K^T 2 ~ S , where K3 and K4 are 
material dependent constants. This model describes the 
change of the exponent of the temperature dependence 
of heat conductivity, as observed in Refs. El IjL Ql , and 
could fit relatively well the experimental data [10|. On 
the other hand, the observation that in narrow bridges 
the exponent of the temperature dependence of the heat 
conductivity apparently converges to 1.5 as the width of 
the bridge decreases [l[ |3j , as well as the increase of the 
ratio k/Cv with temperature Q, cannot be explained in 
this model. 



THE MODEL 

In Refs. p^llll| it was considered simply that the lon- 
gitudinal and transversal polarized phonon modes are 
independently quantized by Neumann boundary condi- 
tions imposed at the membrane surfaces. Nevertheless, 
in a rigorous analysis it must be taken into account that 
the modes with different polarization couple at a free 
surface 0| and because of this, the phonon modes in thin 
membranes are "distorted" and show features that go be- 
yond a simple treatment of the quasi 2D phonon gas. 

For concreteness, we perform calculations for the prac- 
tically important case of SiN^ membranes with thickness 
100 — 200nm having parallel surfaces. Other dimensions 
of the membranes are usually of the order of 100 /urn. 
Consequently, the membranes can be considered as in- 
finitely long and wide. It is further assumed that the 
material is isotropic. The eigenmodes of these kind of 
systems are well known from acoustics and are called 
Lamb waves |l2j . 

Similar to electricity theory, the physical acoustic fields 
can be expressed by a scalar and a vector potential, $ and 
The velocity fields of the longitudinal I and transver- 
sal t phonon modes are then defined as vi = V • $ and 
v t = V xf. Simple wave equations can be derived for 
the potentials, A$ = c~[ 2 d 2 <& and A* = c^ 2 d 2 ^>. Here 
ci and c t {< c;) are the longitudinal and transversal sound 
velocities, respectively. In an infinite three-dimensional 
space this yields one longitudinal and two transversal 
plane wave solutions. 

When however the system is restricted to finite size, 
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FIG. 2: Reflection of elastic waves at the free surfaces of the 
membrane of thickness b. A horizontal shear wave reflects 
also as a horizontal shear wave (a), while a wave longitudi- 
nally or transversally polarized into the xz plane reflects as 
a superposition of longitudinally and transversally polarized 
waves (b). 



the boundaries of the system lead to a coupling of lon- 
gitudinal and transversal modes. This coupling is due 
to the boundary condition at a free surface: the total 
stress should vanish. If the wave incident on the surface 
is polarized along the y direction (see Fig. [5] a), the free 
boundary condition is satisfied if the reflected wave has 
the same polarization. Such a wave is called a horizon- 
tal shear wave (/i-wave). The h- waves do not mix at the 
boundaries with waves of different polarizations and they 
form the typical "box eigenmodes" , with the dispersion 
relation 
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Here fcy is the component of the wave vector parallel to 
the membrane surface and b is the thickness of the mem- 
brane. On the other hand, if the polarization of the in- 
cident wave is either longitudinal or transversal, but in 
the plane xz, then the reflected wave will always be a 
superposition of longitudinally and transversally polar- 
ized waves. These two waves have different propagation 
velocities, and the condition to get eigenmodes is that 
the plane waves reconstruct each other after two reflec- 
tions. The longitudinal and transversal components of 
the eigenmode have the same wave vector component kn 
along the membrane surfaces, but different components 
perpendicular to the surfaces, which we shall call k l ± and 
k^ , respectively. The frequency of the eigenmode is then 
oj 2 = cflk 1 ) 2 = cf(k 1 ) 2 . These equations give a relation 
between the components k l ± and fc^. The eigenmodes 
obtained in this way fall into two classes: symmetric 
(s-wave) and antisymmetric (a- wave), according to the 
symmetry of the velocity field of the wave with respect 
to the plane z = 0. The dispersion relation of the s- and 
a-waves are given by the equations 
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and 



tan(|A;^) _ 

tan(|fci) = '{(k^y-kfy 



(3) 



respectively. The equations (0) and J2J), together with 



<?(*T = c 2 (k 1 ) 2 form a set of transcendental dispers ion relation 



equations, which can be solved only numerically. The 
branches of each dispersion relation are shown in Fig. [3] 
As can be seen, the /i-modes are just usual "box modes" . 
Without mixing at the membrane surfaces, all the modes 
would be of this type and the results would be identical 
to the ones of Refs. [lO^lUl- However, due to the coupling 
of the longitudinal and transversal modes, the dispersion 
relations of the a- and s-modes show some interesting 
properties. 

At fcii = 0, all the excited branches satisfy the relation 
du/dk^ — 0, but for the s- and a-modes, unlike for the 
/i-modes, d 2 u/dk 2 may be either positive or negative. If 
d 2 uj/dk 2 < 0, the dispersion curve has a minimum at 
some value feii > 0. The lowest branch of the /i-modes is 
a straight line, 



Lowest branch of the antisymmetric modes: In the 
case of antisymmetric modes, if (6/2)fcy converges to zero, 
then both (6/2)fc^j_ and [b/2)k l ± converge to zero, but tak- 
ing imaginary values: k\ = ik\ and k l ± = ik l ± , leads to 
a valid solution 0]. Expanding Eq. (0) and using the 
equation ui 2 
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but not for the s- and a-modes. For the s-modes, 
du s $/dk\\ > 0, so the group velocity of long wavelength 
s-phonons is different from zero. On the other hand, for 
the a-modes du} a ,o/dk\\ = and d 2 ui a ,o/dk 2 > 0. There- 
fore the group velocity of the a-modes is zero at long 
wavelength and from this point of view the a-phonons 
are similar to massive particles. 

Since the low temperature thermal properties of the 
membranes are determined by the dispersion relations of 
the lowest branches at long wavelenghts, we shall take a 
closer look at these. 

Lowest branch of the symmetric modes: If (b/2)k\\ 
converges to zero, the solutions (b/2)k t L and {b/2)k l ± of 
Eq. J2J), and satisfying c 2 (fc*) 2 = cf(k 1 ) 2 approach also 
zero, in such a way that k^ is real, while k l ± = ik l ± is 
imaginary [l4|. Then Eq. @ reduces to 
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The solution k^ = is unphysical (does not satisfy the 
boundary conditions if plugged into the stress formulae) , 
so the only solution is ((fc^) 2 — fc|) 2 = 4(fc^) 2 fcjj, which 
yields the dispersion relation 
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So the dispersion relation is linear in the long wavelength 
limit. The group velocity of the s-modes is c s , which, 
since 0.5< 1 — c 2 /cf < 1 |14). takes values between V2c t 
and 2c t . 



(fc') 2 , we get the quadratic 
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as for massive particles of "effective mass 
2c t 6^(c 2 -c 2 )/3c^ 



= h 



(8) 



A plot of the lowest branches of the h-, s- and a-modes 
in the long wavelength limit is shown in Fig. 



THERMODYNAMICS 



Heat capacity 



The qualitative difference between the dispersion re- 
lations of the a-modes, on one hand, and the h- and s- 
modes, on the other hand (see Eqs. 01 El andEJ), have 
important consequences on the thermodynamic proper- 
ties of the membrane. To show this, let us first calculate 
the heat capacity of the membrane. Phonons obey Bose 
statistics, so the average population of each phonon mode 
is n(oj) = 1 / [exp(phuj) — 1]. If we integrate over fcy and 
sum- up the contributions of all the branches, we arrive 
at the expression 
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(exp(/JSw mi « T ) - l) 2 
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where a represents the h-, s-, and a-modes, while Y] m is 
the summation over the branches. The frequencies uj m ,a 
depend also on k\\ ; A is the area of the membrane. 

For uncoupled longitudinal and transversal modes, the 
dispersion laws have the form ui = c a ^J(mr/b) 2 + k 2 . In 

that case the 3D-to-2D crossover in the phonon gas would 
manifest itself through a relatively rapid change of the 
temperature dependence from T 3 to T 2 at the tempera- 
ture T c = Tic t /2bkB, as seen in Fig. |3 (a) 10]. The cor- 
responding asymptotic temperature dependences of the 
heat capacity, following from Eq. 0, are 



m T 3 , T»T C , 
m T 2 , T«T C . 
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FIG. 3: The branches of the dispersion relations of the phonon eigenmodes of a free standing thin membrane, (a) /i-modes, (b) 
s-modes and (c) a-modes 
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FIG. 4: The three lowest branches of the dispersion relations 
of the phonon eigenmodes of a free standing thin membrane 
shown in the limit of long wavelengths. The linear behaviour 
of the h- and s-branches (labeled h and s respectively) as well 
as the quadratic behaviour of the a-branch can be seen here. 

Here, 3/4 = 2/cf + 1/cf, 3/4 = 2/cf + 1/cf, and 
V is the volume of the membrane. The exponent of 
the temperature dependence of the heat capacity, pc = 
<91nCV/91nT, reflects the dimensionality of the phonon 
gas distribution: p c (T > T <: ) = 3 and p c (T < T <: ) = 2 

In the more rigorous case, i. e. when the dispersion 
relations are given by Eqs. JBJ, and J7J, due to the 
quadratic dispersion relation of the lowest o-mode, we get 
a different temperature behaviour for T<T C . Summing 
only over the three lowest branches, we get 

C v » Ak B («T 2 + f3T) 
q= 3^/1 ^ C (2)W 
nh 2 \c 2 c 2 ) nh 2 

(Preliminary estimates show that the contribution of 
dynamical defects to the specific heat are unimportant 
for such thin films; more detailed calculations are in 
progress.) In the high temperature limit, the expected 
T 3 behaviour is obtained, but at temperatures around 
T c , the temperature dependence of pc is somewhat more 



complicated, converging finally to 1, as T — * (see Fig. 
0b). 

Heat conductivity 

A common way to increase thermal insulation of the 
detector mounted on the membrane is to cut the mem- 
brane so that the central part is connected to the bulk 
material only by narrow bridges (see Fig.^l. Therefore, 
another quantity of interest is the thermal conductivity 
(k) along the membrane and bridges. Reported widths 
of the bridges of interest for us are roughly from 4 /im 
upwards 0, Q (We do not discuss here the "dielectric 
quantum wires", like in Ref. EH). As the width of 
the bridge becomes smaller than the phonon mean free 
path in the uncut membrane, the interaction of phonons 
with the bridge edges should become the main scatter- 
ing mechanism. It is very difficult to solve this problem 
either analytically or numerically for the most general 
case. Instead, we shall try to extract the relevant physi- 
cal results for our problem, using a realistic model. The 
cutting process leaves the bridge edges very rough, so we 
shall assume that the phonons scatter diffusively at the 
edges, i.e. scattered phonons are uniformly distributed 
over the angles and branches corresponding to the same 
frequency, ui. 

The general expression for the heat current along the 
rectangular bridge of total length I is 

Q = -\ E M*yM* B )«2(k B )^.g, (ii) 

cr,m,k| 

where r(k||) is the (average) scattering time of a phonon 
having wavevector k|| parallel to the membrane surface 
and belonging to the branch (cr, m). To simplify the writ- 
ing, the dependence on a and m of the quantities in Eq. 
(|llfl was made implicit. The bridge lies along the x direc- 
tion, the z direction is perpendicular to the membrane. 
We denoted by iv,n(k||) = du) a , n /dk\\ the group veloc- 
ity of the phonons. We also assumed that dT/dx is not 
too large, so that the linear approximation, Q oc dT/dx, 
holds. Let us now denote the scattering time of the 
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FIG. 5: The temperature exponents pf = din f /dlnT of the heat capacity (/ = C) and heat conductivity (/ = k) in narrow 
bridges of thin films. The insets show the behaviour of the curves for small values of T/Tc- (a) The temperature exponent of 
CV, as it would be if the modes were not coupled. The dimensionality-crossover around Tc = (TiCi/2/cg) ■ b^ 1 can be seen here 
quite nicely. For a lOOnm thin membrane the critical temperature of SiN x is 237 mK. (b) The temperature exponent of the 
heat capacity of a thin membrane or a narrow bridge, (c) The temperature exponent of the heat conductivity along a narrow 
bridge. 



phonons in the uncut membrane by TM,a,mQz-\\), while the 
scattering time at the bridge edges is te ,a-,m(kii )■ The ef- 
fective scattering time,T - )OT (k|i), is then 



T E,CT,mO £ l| 



Under the assumption of diffusive scattering at the bridge 
ed (k||) can be estimated as r E)0 - im (k||) = 

w / ( v cr,n(k\\) sin??), where z? is the angle between fc|| and 
the cc-direction, and w is the bridge width. If we trans- 
form the summation over k|| in Eq. into an integral 
and write Q = —k, ■ dT / dx, we get 
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where now i M = ^M>,m (k||) = u CT , m (k||) • T M ,a,m(k\\) is 
the mesui free ps-th corresponding to Tyi^^ (k||). Denot- 
ing a = w i 1 Im we can write the integral over $ as 



C{a) = / dti 



COS 2 (#) 

|sin(i9)| + 



= 4 



v 7 1 — x 2 dx 



an 

T 



1 + v 7 ! - a 2 log 



i + vr 



If edge-scattering dominates, i. e. w <C Im, the integral 
only depends logarithmically on a 



C(fc||) « 4 log 



Considering the small bandwidth of occupied phonon 
modes as well as the long mean free paths in the free 
membranes at low temperatures, we can consider C(a) ss 
C = const over the temperature range that we work in. 



With this effective "quasi-constant" mean free path of 
the phonons in the uncut membrane, we can evaluate the 
heat conductivity along the bridge. We did this analyti- 
cally for the low temperature limit of a quasi-2D phonon 
gas and numerically for higher temperatures. The ana- 
lytical low temperature limit gives 
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(13) 



so in the limit of low temperatures k oc T 3 / 2 . The nu- 
merical results are shown in Fig. 0(c). 

If a membrane, which is connected to the bulk ma- 
terial by narrow bridges, is heated by an AC current, 
its thermal cut-off frequency has the expression f c = 
G/Cv = k/(ICv), where k is the thermal conductivity 
of the bridge and Cy is the heat capacity of the mem- 
brane. If all the modes are of the form Q (uncoupling 
of the phonon polarizations at the membrane surfaces), 
then both, Cy and k, ar e propo rtional to T 2 in the limit 
of low temperatures (see |lfj,|ll|) and f c = const. On the 
other hand, using Eqs. (fT77f> and (fC^ we obtain f c oc T 1 ! 2 
in the low temperature limit. The numerical results are 
plotted in Fig. © The increase of f c with T was observed 
experimentally in Ref. 0, but for a wider temperature 
range. 



CONCLUSIONS 

In summary, we used elasticity theory to calculate the 
phonon modes in ultrathin membranes made of homoge- 
nous, isotropic silicon nitride (SiN^). Using the disper- 
sion relations thus obtained, we calculated the heat ca- 
pacity and heat conductivity of the membrane and of 
bridges, cut out of such membranes. In the low tem- 
perature limit the phonon gas becomes two-dimensional 
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FIG. 6: The cut-off frequency of an AC-heated membrane, 
which is connected to the bulk material by narrow bridges. 
Here we normalized the function to its maximum value, so 
that S'c = /c//c,o oc k/(IC v ). For small values of T/T c fc 
becomes linear, which is shown in the inset. 

and populates three branches of the dispersion relations, 
namely the lowest h-, s-, and a-branches (see Fig. [3J. 
At low temperatures, the dispersion relation correspond- 
ing to the lowest a-branch is quadratic in fen , while the 
other two are linear, so at low temperatures the a-branch 
gives the dominant contribution. Quite surprisingly, this 
implies that the universal behavior of heat capacity in 
two-dimensional systems is obeyed also by the phonon 
gas at low temperatures, where CV °c T (see Ref. [lif 
and citations therein). 

In the calculation of thermal conductivity along 
bridges, we assumed diffusive scattering of phonons at 
the bridge edges. This is justified by the fact that the cut- 
ting process leaves these edges very rough. If the width 
of the bridge decreases below a certain value (depending 
on the mean free path of the phonons in the uncut mem- 
brane) the interaction with the edges becomes the main 
scattering process for the phonons, see Eq. I|12|) . In such 
a case, at low temperatures it was found that k oc T 3 / 2 
(Eq.EU. 

If a membrane, connected to the bulk material by nar- 
row bridges, is heated by an AC current, the amplitude of 
the temperature oscillations in the membrane has a cut- 
off around the frequency f c = G jCy . In the low temper- 
ature limit, this cut-off frequency shows an increase with 
the temperature, as T 1 / 2 . Preliminary experimenatal re- 
sults show an increase of f c with T, but on a temperature 
range much wider than the one in Fig. |SJ This seems to 
suggest that at higher temperatures other processes have 



to be taken into account in the calculation of thermal 
characteristics of SiN x membranes. 

ACKNOWLEDGEMENTS 

We thank Prof. J. Lothe for various discussions about 
elasticity theory and Dr. I. M. Maasilta for useful com- 
ments on the manuscript. We also acknowledge the fi- 
nancial support from the Academy of Finland under 
the Finnish Center of Excellence Program 2000-2005 
(Project No. 44875, Nuclear and Condensed Matter 
Physics Program at JYFL) and the NorFA Networks, 
"Mesoscopic Solid State and Molecular Electronics" and 
"Quantum Transport in Nanoscale systems" . 



[1] M. M. Leivo and J. P. Pekola, Appl. Phys. Lett. 72 1305 
(1998) 

[2] A. Luukanen, M. M. Leivo, J. K. Suoknuuti, A. J. Man- 
ninen, and J. P. Pekola, J. Low. Temp. Phys. 120, 281 
(2000). 

[3] J. P. Pekola, A. J. Manninen, M. M. Leivo, K. Aru- 

tyunov, J. K. Suoknuuti, T. I. Suppula and B. Collaudin, 

Physica B 280, 485 (2000). 
[4] T. S. Tighe, J. M. Worlock and M. L. Roukes, Appl.Phys. 

Lett. 70 2687 (1997) 
[5] W. Holmes, J. M. Gildemeister, P. L. Richards, and V. 

Kotsubo, Appl. Phys. Lett. 72, 2250 (1998) 
[6] A. L. Woodcraft, et al., Physica B 284-288, 1968 (2000). 
[7] S. Hunklinger and A. A. Raychaudhuri, in Progress in 

Low Temperature Physics, vol. IX, ed. by D. F. Brewer 

(Amsterdam, Elviser), 267 (1986). 
[8] S. Wenzel, Ph.D. thesis, University of California at 

Berkeley, 1992. 

[9] J. M. Ziman, Electrons and Phononos: The Theory of 
Transport Phenomena in Solids, Oxford at the Clarendon 
Press, 1960. 

[10] D. V. Anghel, J. P. Pekola, M. M. Leivo, J. K. Suoknuuti 

and M. Manninen, Phys. Rev. Lett. 81, 2958 (1998) 
[11] D. V. Anghel, Phys. Rev. B 59, 9854 (1999) 
[12] B. A. Auld, Acoustic Fields and Waves in Solids, volume 
II, second edition, Robert E. Krieger publishing company 
1990. 

[13] L. D. Landau and E. M. Lifshitz, Theory of Elasticity, 

Pergamon Press, 1959. 
[14] T. H. Kiihn, phonons in nanoscopic films, masters thesis 

2002. 

[15] L. G. C. Rego and G. Kirczenow, Phys. Rev. Lett. 81, 232 
(1998); D. E. Angelescu, M. C. Cross, and M. L. Roukes, 
Superlattices and Microstructures 23, 673 (1998); K. 
Schwab, E. A. Henriksen, J. M. Worlock, and M. L. 
Roukes, Nature 404, 974 (2000). 

[16] D. V. Anghel, J. P hys. A: Math. Gen. 35, 7255 (2002); 
cond-mat/0105089 



